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Article history: The mean illumination conditions and surface temperatures over one orbital period are calculated for the Asteroid 4 

Received 26 July 201 1 Vesta using a coarse digital elevation model produced from Hubble Space Telescope images. Even with the anticipated 
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Accepted 4 November 2011 currently exist on Vesta due to its large axial tilt («27°). However, under present day conditions, it is predicted that 
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modeling of main belt asteroids, suggests that water ice could survive in the top few meters of the vestal regolith 
on billion-year timescales. 
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1. Introduction 

As with other airless bodies in the Solar System, the surface of Vesta is directly 
exposed to the full solar spectrum. The degree of solar illumination plays a major 
role in processes at the surface, including heating, electric charging, chemistry, 
and exospheric production. It is likely that illumination (surface temperature) con- 
trols the transport and deposition of volatiles at Vesta, as has been proposed at the 
Moon (e.g., see Sunshine et al., 2009; Paige et al., 2010). If there are permanently 
shadowed regions (PSRs) on Vesta - or at least regions that, on average, receive very 
little solar illumination - then they could act as cold traps for volatiles. Schorghofer 
(2008) estimates that water ice on main belt asteroids can persist within the top 
few meters of a dusty surface layer over billions of years, if the mean surface tem- 
perature is ;$145 K. 

We predict the mean illumination conditions and surface temperatures over a 
vestal year (1325 days) using a modified version of the Vesta digital elevation mod- 
el (DEM) produced by Thomas et al. (1997b). This allows us to make an assessment 
of the ability of volatiles, in particular water ice, to survive in the vestal regolith 
under present day conditions. The general methodology employed here with a 
relatively coarse DEM could readily be applied to other airless bodies in the Solar 
System. Determining illumination and thermal conditions is particularly timely 
given the Dawn mission’s arrival at Vesta (Russell et al., 2007), with a payload 
including the Gamma Ray and Neutron Detector (GRaND). The regions with mean 
surface temperatures ;$145 K predicted here could correlate with possible future 
detections of hydrogen-rich deposits by GRaND, which may be indicative of 
sub-surface water ice. 
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2. The orbital and shape characteristics of 4 Vesta 

Vesta is located in the inner main asteroid belt, and its current orbital charac- 
teristics are summarized in Table 1. Based on Hubble Space Telescope (HST) images 
acquired in 1994 and 1996, Vesta is known to be a triaxial ellipsoid with radii of 
289, 280 and 229 ±5 km with a spin pole solution of RA = 301°, Dec = 41°, ±5° 
(revised values from Thomas et al. (1997b)). This corresponds to a large axial tilt, 
or obliquity, of 27.2 ± 5°, which is in good agreement with more recent studies (Li 
et al., 201 1 ). Vesta has a prograde rotation with a relatively short spin period for 
an asteroid of 5.342 h. Vesta is close to being an oblate spheroid, but for the 
enormous impact crater centered near its south pole with a diameter of ~460 km 
(Thomas et al., 1997a). 

For the shape of Vesta we use the 5x5° DEM based on HST images (Thomas 
et al., 1997b). The orbital information for Vesta is obtained from the Dawn SPICE 
kernel, which is currently limited to the 1900-2100 epoch. However, before this 
DEM could be used for predicting illumination conditions, some modifications 
were required in order to produce more realistic results. Firstly, we had to smooth 
out some irregularities in the DEM poleward of 80°N/S that were clearly associ- 
ated with the meridional convergence, as opposed to any real topographic fea- 
tures. We found that this could be effectively dealt with by simply smoothing 
the relevant data points on the original DEM, which was achieved by calculating 
an average radial distance using one data point either side in latitude (±5°) and 
four points either side in longitude (±20°). Secondly, the 5x5° DEM is too coarse 
for producing realistic illumination predictions («22.5 km at the equator): so we 
decided to interpolate the DEM to a 1 x 1° resolution («4.5 km) using the modi- 
fied Shepard’s method (Renka, 1988), which is a multivariate interpolation tech- 
nique that employs an inverse distance weighting algorithm. To our knowledge, 
this is the best interpolation method for increasing the resolution of coarse DEMs, 
since it has the advantages of producing smooth results while still exactly match- 
ing the values of the original DEM. Using this modified DEM, we produced revised 
estimates for the mean radius, surface area, volume and mean density of Vesta, as 
detailed in Table 1. 
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Table 1 

Summary of the orbital, physical and thermal properties relevant to Vesta. 


Property 

Value 

Aphelion 3 

2.572 AU 

Perihelion 3 

2.151 AU 

Eccentricity 3 

0.08890 

Semi-major axis, 3 D 

2.36158 AU 

Inclination to ecliptic plane 3 

7.13404° 

Orbital period 3 

1325 days (3.63 years) 

Precession period 3 

81,730 years 

Spin period 3,13 

5.342 h 

Spin pole solution 13 

RA = 301°, Dec = 41°, ±5° 

Obliquity (axial tilt) b - c 

27.2 ± 5° 

Triaxial ellipsoidal radii b 

289, 280, 229 ± 5 km 

Mean radius d 

266.6 ±8.8 km 

Surface area d 

8.79 ± 0.26 x 10 5 km 2 

Volume d 

7.81 ± 0.26 x 10 7 km 3 

Mass e 

2.586 ± 0.002 x 10 2o kg 

Average density d e 

3.31 ±0.11 x 10 3 kgnr 3 

Solar constant, S 0 

1366 Wm- 2 

Solar constant at Vesta, S 0 /D 2 

245 W irr 2 

Bond albedo/ A 

0.20 

Thermal emissivity/ e 

0.9 

Diurnal skin depth/ S d 

0.00192 m 

Seasonal skin depth/ <5 S 

0.148 m 

Mean nightside surface temperature (hybrid 
model) 

80 K 


a Bills and Nimmo (2011). 
b Thomas et al. (1997b) revised values. 
c Li et al. (2011). 

d The modified Vesta DEM presented here. 
e Kuzmanoski et al. (2010). 
f Chamberlain et al. (2007) and references therein. 


3. Illumination model 

The illumination model used here has been optimized for speed, accuracy and 
versatility by combining the ray tracing approach (e.g., Bussey et al., 2010) with 
the “horizon method” described by Mazarico et al. (2011). Since we are just inter- 
ested in the illumination conditions for this study, we only need to calculate: (1 ) the 
fraction of the solar disk visible from the surface (accounting for topography along 
the horizon); (2) the solar incidence angle (including the effect of local slope); and 
(3) the normalized solar flux, which is determined using the previous two parame- 
ters. Given the uncertainties in the DEM, we do not consider sunlight scattered from 
other surface elements, just direct solar illumination. 

Similar to Mazarico et al. (201 1 ), the first step is to generate a database contain- 
ing the elevation of the horizon at 36 equally-spaced azimuthal directions for every 
point on the DEM. Using this database, the elevation and slope of the horizon can 
then be calculated very efficiently for any given azimuthal direction. The elevation 
and angular extent of the Sun is given by the SPICE kernels, which then allows the 
fraction of the solar disk visible from the surface to be calculated. The DEM surface 
normals are calculated by taking the cross product of adjacent grid points, which 
then permits the solar incidence angle to be determined. However, unlike Mazarico 
et al. (2011), we do not rely on gnomonic projection in our calculations; therefore, 
we are not restricted to the polar regions and are able to achieve global coverage. 
Our model has been validated by comparisons with illumination simulations of 
the lunar polar regions generated by Mazarico et al. (201 1 ) using Lunar Orbiting La- 
ser Altimeter (LOLA) DEMs. 

In our simulations we consider the average illumination conditions over one 
vestal year, which should be indicative of conditions during at least the past few 
hundred thousand years (cf. Bills and Nimmo, 201 1 ). We found that the normalized 
solar flux results were typically accurate to within <0.01% when calculations were 
performed 72 times per rotation period at 360 locations along the orbital path (sep- 
arated equally in time). 

4. Hybrid thermal model 

In order to use our illumination predictions to assess the ability of volatiles to 
survive in the vestal regolith, we need to employ a thermal model. Schorghofer 
(2008) made theoretical estimates of the loss rate of buried water ice from main 
belt asteroids. He proposed that water ice can survive within the top few meters 
of a dusty surface layer for billions of years, if the mean surface temperature (7) 
is <145 ± 10 K. As part of his study, two end member thermal models were investi- 
gated: the fast-rotating or isothermal latitude model, which assumes infinite ther- 
mal inertia; and the non-rotating or standard thermal model, which assumes zero 
thermal inertia (see Schorghofer (2008), and references therein). When compared 


with a model that included a more realistic multi-layered treatment of thermal 
inertia, it was shown that neither of the end member models provided reasonable 
estimates for (T). Therefore, we employ a hybrid thermal model that uses the illu- 
mination predictions to estimate the mean dayside temperatures ( T d ), and assumes 
a constant mean nightside temperature (T n ). The assumptions that the dayside is in 
thermal equilibrium with the incident solar flux, while the nightside temperature 
does not vary significantly, is consistent with thermal observations of other airless 
bodies (e.g., Paige et al., 2010). The mean dayside temperature for each surface ele- 
ment is determined using the Stefan-Boltzmann law, (T d ) 4 = (1-A)(S 0 (F)/£<jD 2 (V)); 
where o is the Stefan-Boltzmann constant, (V) is the average fraction of the solar 
disk visible from the surface, (F) is the average normalized solar flux scaled to the 
semi-major axis D (see Table 1 for more details). The overall mean surface temper- 
ature for the hybrid model is then given by (Th) = {V)(T d ) + { 1 - ( V))(T n ), which 
produces estimates that lie between the two end member models, as expected. 

To verify this hybrid approach we compare with results from Schorghofer 
(2008) that indicated «20% of the surface area of a spherical body with Vesta’s 
semi-major axis and obliquity would have (T) < 145 K, which matches our model 
if (T n ) = 83 K. In Schorghofer’s calculations, the orbital and thermal characteristics 
used are broadly consistent with those expected for Vesta (e.g., see Table 1), apart 
from a lower Bond albedo of A = 0.05 (used in the above comparison). In our study 
we use A = 0.2 and ( T n ) = 80 K with the hybrid thermal model, which for a Vesta-like 
spherical body results in (T/,) < 145 K over %40% of its surface. We estimate the (T/,) 
uncertainties to be typically <10 K. For context, these thermal characteristics for 
Vesta predict a subsolar surface temperature at the semi-major axis of 249 K. 


5. Results 

Illumination and thermal modeling predictions are presented in planetocentric 
coordinates with west longitudes plotted to the right of the prime meridian in order 
to be consistent with earlier presentations of Vesta DEMs (Thomas et al., 1997a,b). 

5.1. Average illumination conditions over a vestal year 

We present global predictions for the average illumination conditions at Vesta 
during one vestal year, with particular attention paid to the southern polar region 
(Fig. 1). The top row shows the modified Vesta DEM (1 x 1° resolution), which 
clearly reveals the large topographic variation relative to Vesta’s basic triaxial ellip- 
soid shape. The giant impact crater can be readily identified in the southern polar 
region with its central peak at 65°S, 104°W and the lowest point at 71°S, 111°E. 
In contrast, the northern polar region appears to be relatively featureless. 

The second row shows the mean solar visibility (V), which reveals that the solar 
disk is typically seen about half the time in the tropics (between the dotted lines), as 
expected for a diurnal variation. At mid- to far-northern latitudes, ( V) tends to in- 
crease slightly, probably due to Vesta’s irregular oblate shape. Conversely, in the 
southern hemisphere, the effects of the giant impact crater are clearly evident with 
many distinct regions of low (V). As expected, the point with the highest (V) (0.566) 
is located close to the central peak at 65°S, 100°W; while the point with the lowest 
(V) (0.363) is close to the topographic low at 70°S, 108°E. The third row shows the 
mean solar incidence angles (/), where the expected general trend of low (I) around 
the equator, which then increases toward the poles, is evident. Due to the promi- 
nence of the central peak, the (I) values are typically much lower there than for 
the rest of the southern polar region. 

The maps of mean normalized solar flux (F) (fourth row) have been scaled to the 
semi-major axis, and so can be used directly with the hybrid thermal model. The 
highest (F) value (0.302) is predicted at the equator at 0°, 4°E; this is consistent with 
expectations since each surface element will be illuminated about half the time, and 
so combined with a simple cosine dependence for a spherical body gives a maxi- 
mum value of 0.318 (l/7t). The lowest (F) value (0.145) is located within the giant 
crater at 73°S, 30°E. In spite of this, the northern and southern hemispheres have 
roughly the same surface area receiving (F) < 0.15. As anticipated for a body with 
such a significant axial tilt, we do not identify any permanently shadowed or per- 
manently illuminated regions on Vesta. However, we have not yet considered the 
effects of finer-scale topography (see Section 5.2). 

Using the hybrid thermal model described in Section 4 we convert the mean 
illumination conditions to mean surface temperatures (T/,) (bottom row), which 
gives a range of values from 123.4 to 150.8 K, with a global annual average of 
142.7 K. Based on our results, under present day conditions about 52% of Vesta’s 
surface is cold enough for water ice to have persisted in the top few meters of 
the regolith for billions of years should it have been delivered there. The 
(T/,) = 145 K contour is indicated by the gray solid lines. 

5.2. Influence of finer-scale topographic features on illumination and thermal conditions 

Since the original DEM used in this study has a resolution of 5 x 5°, the modi- 
fied version is still only able to really resolve topographic features with spatial 
scales of ~100km. It is well known that topography with much smaller spatial 
scales than this can have a significant influence on illumination conditions (e.g., 
Mazarico et al., 2011). Therefore, to investigate this influence we convolved a 
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Fig. 1. Predictions for the illumination conditions and surface temperatures at Vesta, averaged over one vestal year, based on a modified version of the Thomas et al. (1997b) 
DEM. Shown are global predictions in equidistant cylindrical projection (left column), and predictions for the southern polar region in orthographic projection (right column). 
Also indicated are the equator (black dashed line), tropical circles (dotted black lines at 27.2°N/S), and polar circles (white dashed circles at 62.8°N/S). Shown from top row: 
the radial planetocentric distance of the surface; the visible fraction of the solar disk; solar incidence angle; normalized solar flux (where 1 corresponds to the Sun directly 
overhead at the semi-major axis); mean surface temperature derived using our hybrid thermal model. (For interpretation of the references to color in this figure legend, the 
reader is referred to the web version of this article.) 
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Fig. 2. Predictions for the influence of finer-scale topography on the illumination conditions and surface temperatures at 4 Vesta, averaged over one vestal year. Shown are 
predictions in orthographic projection for simple craters with diameters of 1 km (left column) and 10 km (right column) centered at 73°S, 30°E. Parameters shown are the 
same as in Fig. 1. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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simple crater shape with the modified DEM centered at 73°S, 30°E, the location pre- 
dicted to receive the minimum (F) in Section 5.1. In this study we use the typical 
shape characteristics for simple lunar craters with diameters <15 km, as described 
by Horz et al. (1991) and references therein. In order to assess the dependence on 
the spatial scale of the craters, we investigate the illumination conditions and sur- 
face temperatures for a 1 km and a 1 0 km diameter crater averaged over one vestal 
year. 

The results for the 1 km diameter crater are shown in the left column of Fig. 2. 
Despite only producing an apparently minor perturbation to the local topography, 
this crater has a significant effect on average illumination conditions and surface 
temperatures. The values of (F) at the center of the crater floor drop by about a fac- 
tor of 2 to ^0.08, which results in (T h ) values of only fslOl K. This is almost 30 K 
colder than if there were no crater present, and could result in water being ther- 
mally stable on the surface at that location for much of the vestal year (cf., Paige 
et al., 2010). The results for the 10 km diameter crater, shown in the right column 
of Fig. 2, reveal that this size crater produces a significant perturbation to the sur- 
rounding topography. In this example we can see how (V) values at the center of the 
crater decrease to ^0.13 (almost a factor of three drop), but increase along the nar- 
row edge of the crater rim to ~0.5. Similar to the 1 km crater case, there is a signif- 
icant decrease in (I) values on the equatorward facing southern crater wall of up to 
«15°. This suggests that scattered sunlight could sometimes play an important role 
in the thermal conditions within such a crater (cf. Paige et al., 2010). The values for 
(F) decrease to ?s0.065 on the crater floor (more than a factor of two drop), while 
they increase to ^0.18 along a narrow crescent on the southern wall (up by 
«20%) due to the decrease in (I) values discussed previously. The predicted change 
in (Th) ranges from a decrease to ^97 K on the crater floor, to an increase along the 
edge of the crater rim to «130 K. Despite the increased (F), the southern crater wall 
is colder on average due to the decrease in (V). Conversely, the sharp increase in (T h ) 
along the edge of the crater rim is due to the increased (V). 

Neither a 1 km nor a 10 km diameter crater is able to produce a PSR at this loca- 
tion, although the latter gets very close. In both cases there is a large perturbation to 
the average illumination conditions that produces an overall decrease in (Th), which 
is anticipated to significantly affect the thermal stability of volatiles at the surface, 
as well as those buried in the regolith. 

6. Discussion and conclusions 

The depth to which the solar thermal perturbation can penetrate into the reg- 
olith can be assessed using the thermal skin depth <5, which is dependent on the 
time period of the temperature oscillation (see Chamberlain et al. (2007), and ref- 
erences therein). Since Vesta has an obliquity of «27.2°, we must consider both 
the diurnal skin depth S d (rotation period), and the seasonal skin depth <5 S (orbital 
period). The S d value is relevant to equatorial and mid-latitude regions; while S s 
is applicable to the polar regions, which can be exposed to near-constant solar illu- 
mination during their summer seasons (see Table 1). If water ice does exist in the 
vestal regolith, then it will likely be covered by a desiccated layer of at least a 
few S in thickness (cf. Chamberlain et al., 2007). Therefore, in the polar regions of 
Vesta it seems reasonable to expect any significant abundance of water ice to be 
buried at least a few tens of cm below the surface. 

Recent observations of a 3 pm absorption feature on the Moon with band 
depths of up to ^15% indicate that its entire surface is hydrated during some por- 
tion of a lunation, with diurnal variations and other characteristics possibly consis- 
tent with hydroxyl formed by the interaction between solar wind protons and 
oxygen-bearing minerals (e.g., Sunshine et al., 2009). The 3 pm absorption feature 
has also been observed on Vesta (Hasegawa et al., 2003), but the band depth was 


only «1% which seemed too weak to be consistent with the above solar wind 
implantation hypothesis. However, the apparent lack of space weathering on Vesta 
suggests that it is sufficiently magnetized to be able to shield most of its surface 
from solar wind protons (e.g., Vernazza et al., 2006). This proposed proton-shielding 
could also be the reason for the relatively weak 3 pm absorption feature, and be 
indicative of lower than expected abundances of hydrated and/or hydroxylated 
minerals on the surface of Vesta. 

Based on our illumination and thermal predictions for present day conditions, it 
is possible that water ice could have survived in the top few meters of the regolith 
in the polar regions of Vesta for billions of years (over about half of Vesta’s surface) 
should it have been delivered there. As expected, topography plays an important 
role in determining the average illumination conditions and surface temperatures. 
In particular, small scale topographic features, such as ~1-10 km diameter craters, 
can have a significant effect, and perhaps even permit water ice to be thermally sta- 
ble at the surface ( < 100 K) in some polar locations for much of the vestal year. 
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